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Abstract. - It is shown that the elimination of the discrete transverse motion in a waveguide 
of arbitrary shape may be described in terms of a non-abelian gauge field for the longitudinal 
dynamics. This allows for an exact treatment of the scattering between different modes by elim- 
inating the gauge field at the expense of a non-diagonal matrix of local subband energies. The 
method is applied to calculate the local density of states (LDOS) in a quantum point contact. 
Contrary to the total conductance which is well described by an adiabatic approximation, mode 
mixing turns out to play a crucial role for local properties like the LDOS. 



The discovery of quantized longitudinal conductances in quantum point contacts created 
in a two-dimensional ballistic electron gas || has launched a great number of theoretical 
papers dealing with quantum mechanical scattering in waveguides with a general wide-narrow- 
wide type geometry ||. Precise numerical calculations of the transmission amplitudes t mn 
between the transverse modes m and n in such a structure have shown jij that the adiabatic 
approximation t mn ~ S mn ~ while not exact - works extremely well for the linear conductance 



of the point contact at zero temperature. The effect of scattering between different transverse 
modes is therefore negligible for the conductance, thus justifying the adiabatic approximations 
used in the early explanations ^ |(| of the observed quantization of G. In our present work, we 
develop a novel method dealing with quantum mechanical scattering in waveguides of arbitrary 
width. It is based on a gauge theoretic description of the effect of transverse modes on the 
dynamics of the longitudinal motion, which allows for an exact treatment of intersubband 
scattering. Using this formalism and the recursive Green function method Q, we calculate 
the local density of states in a quantum point contact. In the vicinity of the constriction this 
local property is strongly affected by the scattering between different transverse modes. 

We consider a general wave guide in two dimensions with an arbitrary but fixed confin- 
ing potential V(x,y). To determine the exact scattering wave functions of the stationary 
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Schrodinger equation 



2M 



V 2 + V(x,y) 



ip(x,y) = Eip(x,y) 



(2) 



we first split off the discrete motion in the transverse coordinate y. Quite generally, for 
any confining potential V(x,y) which approaches infinity as y — > ±00, the one-dimensional 
Schrodinger equation for the transverse motion has a complete and discrete basis of eigenstates 
&n(y) f° r every value of the longitudinal coordinate x. These eigenstates and eigenenergies 
are of course position dependent. The position dependent eigenenergies e n (x) give rise to 
potential barriers for the electrons traversing the wave guide structure in longitudinal direction. 
The corresponding eigenfunctions can be chosen real globally, which is always possible in the 
absence of a magnetic field. The completeness of the transverse basis allows one to expand 
the total wave function in the form if)(x,y) = ^ n Xn( x )^n(v)- Inserting this ansatz into the 
two-dimensional Schrodinger equation, we obtain 



X(x) = Ex(x). 



(3) 



Here x{x) = (\i(x), Xi{ x )i ■ ■ ■) IS the vector of the longitudinal wave functions for the different 
transverse eigenstates n and p — —ihd x is the momentum operator for the longitudinal motion. 
To simplify the notation we have introduced the following matrices 



+00 

A nm {x) := l dy ®Z(y) 



dx 



and B r , 



+00 

t {x) := / dy $l{y) 



dx 2 



(4) 



and also the diagonal matrix of the local subband energies 

2M dy 2 



i(x) := (n{x)\ 



V{x,y)\m(x)) = e„(x) ■ S n 



(5) 



It is straightforward to see that A is antisymmetric, since the wave functions are real valued 
everywhere. Obviously A(x) and B{x) are smooth for confining potentials, which are smooth 
functions of the longitudinal coordinate x. If there are regions in the form of an ideal lead, we 
have A{x) = B{x) — and the above system decouples trivially. 

A simple calculation leads to the following relation between the matrices A(x) and B(x) 



B(x) 



dA(x) 



dx 



+ A 2 {x), 



which was first discovered by Kuperin et al. 
the form 



. This relation allows one to rewrite eq. 



1 

2M 



(p-ihA(x)) + e(x) 



X(x) = Ex(x). 



(6) 



(7) 



Formally this is an tensorial Hamiltonoperator acting on an infinite-dimensional vector of wave 
functions. Apparently, the form of our Hamiltonian now resembles the quantum mechanical 
Hamiltonian of a charged particle in a magnetic field, however, in the present case the gauge 
field A is a tensor instead of a vector. 

We now consider unitary transformations performed on the set of transverse eigenfunctions, 
i.e. we introduce a local gauge transformation (the bar denotes complex conjugation) 



\k(x)) 



S k n{x)\n(x)), 



(8) 
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which mixes the different transverse eigenmodes. Since we do not want to change the physics 
via this operation, we also have to change our longitudinal coefficients according to 



(9) 



leaving the entire wave function invariant. It is straightforward to show, how the differential 
equation ^ for the longitudinal coefficients x(x) changes under this transformation. The 
transformation law for the diagonal matrix of subband energies e(x) is simply 



i(x) = S(x)e(x)S + (x). 



(10) 



The new and old gauge potentials A(x) and A{x), however, are related by the more complicated 
transformation law 



A{x) = S(x)A{x)S + (x) - ^^S+(x), 

ax 



(11) 



well-known from non-abelian gauge theories The transformed Schrodinger-equation for 
the longitudinal coefficients then has the same form as (R), i. e. 



1 / - \ 2 

— [p - ihA{x) j + g(x) 



X(x) = Ex(x). 



(12) 



Since i}}(x.y) is unchanged, all observables are invariant under this gauge transformation. 
Choosing a gauge in this context means to choose a certain basis of transverse wave functions 
at every space point x. 

The formalism above can now be used to describe the mixing between the different transverse 
modes - which is formally described by the gauge potential A(x) - in a very explicit form, that 
also turns out to be convenient for a numerical treatment. To this end a gauge transformation 
is performed to a gauge in which A(x) vanishes identically. From ([ll]) this requires that 
S(x)A(x) — dS(x)/dx or 



dS+(x) 
dx 



S(x) = -A(x). 



This differential equation for the transformation matrix S(x) can be solved formally by 



S(x, xq) = R 



exp / d£ A(0 



(13) 



(14) 



where R is the space-ordering operator, defined in complete analogy to the time-ordering 
operator known from time-dependent perturbation theory. The antisymmetry of the gauge 
potential A ensures the unitarity of the evolution operator S(x,Xq). In order to understand 
the physics behind our special gauge, we note that A is defined by 



d_ 

dx 



n(x)) 



y^ j A nm (x)\m(x)). 



(15) 



The transverse basis in the gauge where A(x) = is therefore obviously the one with a fixed 
reference point xq. Thus by expanding the total wavefunction in terms of transverse modes 
which are independent of the longitudinal coordinate x, the gauge potential clearly vanishes. 
The fixed basis set $^°(y) is related to the locally adiabatic basis $n(v) & 1 anv point x by the 
unitary transformation S(x, xq) as given in (fbij). As a result, the original diagonal matrix s(x) 
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of the adiabatic subband energies is transformed to a non-diagonal form via the transformation 
law (|l0|). Mode-mixing is thus no longer described by a non- vanishing gauge potential but is 
instead due to the off-diagonal terms in the transformed energy matrix e(x) . 

For an application of the general theory above, it is necessary to specify the confining 
potential. In order to obtain analytical results as far as possible, we choose a harmonic 
confinement of the form 

n*.v) = ^(tMV- (16) 



2 u \b(x 

The effective width b(x) has a minimum bo = b(x = 0) at the origin and widens to a finite value 
boo = b(x — > ±oo) at \x\ larger than a characteristic constriction length L. Quite generally, for 
a harmonic confining potential, the gauge potential A(x) may be written in the form 

("0 

where a is the usual harmonic oscillator destruction operator and the prime denotes differenti- 
ation with respect to x. Since A(x) factorizes into a simple function of x and a fixed operator, 
the spatial ordering symbol in ( |l4| ) is irrelevant in the present case. The evolution operator S 
is therefore obtained in explicit form as 



S(x, Xq) = exp 



1 fb(x) 



4 \b(x ) 



(a +2 -a 2 ) 



(18) 



This operator turns out to be identical with the so called squeeze-operator well known in 
quantum optics || . Its matrix elements in the transverse oscillator eigenstates can be expressed 



analytically 1 1 in terms of hypergeometric functions, which are not given here explicitly, 
however. 

In order to solve the coupled differential equations for the longitudinal wave functions x(x), 
we eliminate the gauge potential and discretize the transformed equation in steps of size a, 
such that x( x = * a ) = Xi- Transforming back to the original adiabatic basis with a diagonal 
matrix e{x) of subband energies, the resulting equations are 

g t + 2 J Xi ~ Si+i^Xi+i - Si-i,iXi-i = Exi (19) 

in an obvious notation for e ., Si-fi,j and with dimcnsionless energies in units of Ti 2 /{2Ma 2 ). 
This is a tight-binding like model in one dimension with as many orbitals per site i as number 
of transverse modes included (about 20 in the numerical calculations). A convenient and 
numerically stable method to solve it is the recursive Greenfunction technique Q . This allows 
one to calculate the matrix of the local Greenfunctions G mn (i,j) at a given energy ep from 
which the local density of states 

p(x, y;e F ) = --Y] $^(y)$^(y) • Im G mn {x, x) (20) 

mn 

is finally obtained. This quantity can be directly measured in an experiment with a scanning 
tunneling microscope (STM), provided the two-dimensional electron gas can be realized on a 
free surface like in InAs-systems. Indeed the local tunneling current between an STM tip and 
a conducting sample is directly proportional to p(x\£f) ]TT[ . In the standard case where the 
electron gas resides below a surface layer, the LDOS may be determined via a local capacitance 
spectroscopy as was recently demonstrated by Ashoori [Q . An example for the local density 
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Fig. 1. - Local density of states for a quantum point contact with a harmonic confinement of width 
b(x) as indicated. In the vicinity of the constriction the "exact" calculation (left) differs strongly from 
the result obtained in an adiabatic approximation (right). 

of states is shown in fig. 1 for a smooth boundary function b(x) with a width ratio bo /fooo = 0.1 
as indicated in the figure. The characteristic length L is 250 nm while the Fermi energy 
£p = Hloq/2 is choosen such that it coincides with the lowest transverse mode at the narrowest 
point x — 0. Including the quantum mechanical tunneling and reflection over the barrier 
the associated linear conductance is equal to e /h, i.e. we are at the inflection point of the 
transition between pinch off and the first conductance plateau at G = 2e 2 /h. This situation is 
nicely reflected in the exact calculation of the LDOS shown in fig. 1 on the left, where a narrow 
dip at the constriction is just about to be closed. In the wide regime where four transverse 
modes are occupied, the LDOS exhibits a complex interference pattern reflecting the square 
of the exact wavefunction ^2 n Xn(x)^„(y) at the Fermi energy. In this asymptotic regime the 
structure is reproduced rather well in an adiabatic approximation where all matrices S in ( |l9| ) 
are replaced by a unit matrix. The resulting LDOS is shown in fig. 1 on the right for the same 
parameter values. It is very similar to the exact result in the regime where b(x) is constant. 
Near the constriction, however, the adiabatic approximation obviously differs strongly from 
the exact result and in particular it fails to reproduce the almost pinch-off structure. Mode 
mixing therefore plays a crucial role for local properties like the LDOS. 

In conclusion we have developed a novel method to treat quantum scattering in general 
wave guides via a gauge field approach and have applied it for a calculation of the LDOS in 
a quantum point contact. The nontrivial local structure found here may possibly be observed 
using modern scanning probe techniques. Our method can also be used to determine the local 
potential distribution in the presence of a finite current [|l3j and the local emissivities and 
injectivities which determine the ac- transport properties at low frequencies Jl4[ . 

It is a pleasure to acknowledge useful discussions with Bert Lorenz and David Wharam 
about the possibilities of measuring the LDOS in a quantum point contact. This work was 
supported by the SFB 348 Nanometer-Halbleiter-Bauelemente. 
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